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V^'. Abstract 

' The STAR collaboration at RHIC is measuring the production of electron- 

^SJ ■ positron pairs at small impact parameters, larger than but already close to the 

i range, where the ions interact strongly with each other. We calculate the total 

^ cross section, as well as, differential distributions of the pair production process 

with the electromagnetic excitation of both ions in a semiclassical approach and 
, within a lowest order QED calculation. We compare the distribution of electron and 

positron with the one coming from the cross section calculation without restriction 
on impact parameter. Finally we give an outlook of possible results at the LHC. 

> : 

O ■ 1 Introduction 
(N 

^ I Pair production in relativistic heavy ion collisions has attracted interest in the past mainly 
O ' due to the fact that the strong fields allow for multiple pairs to be produced. At impact 
parameter of the order of twice the nuclear radius, but still larger than this, so that the 
two ions do not interact hadronically with each other (that is the regime of the so-called 
O I "ultraperipheral collisions" UPC), the total pair production multiplicity is found to be 

about 1.5 for AuAu collisions at RHIC and about 3.9 for PbPb collisions at the LHC. 
These results are based on a lowest order QED calculation [T] and one might suspect that 
at these small impact parameters the strong fields of the two ions do lead to higher order 
corrections. 

^ I Therefore it is of interest to measure electron-positron pairs produced in such collisions 

and compare their distribution with theoretical predictions, e.g., in lowest order QED. The 
STAR collaboration has recently measured e~^e~ pairs in collisions, which were selected by 
a trigger, looking for the simultaneous excitation of the two ions (mainly to the GDR) in 
addition to the pair production process [3 El E! , see Fig. [H Such an event is characterized 
by the subsequent emission of one or only a few neutrons, which are then detected in the 
forward ZDC ("Zero Degree Calorimeter"). This tagging was first proposed in order to 
study vector meson production E] • Mutual excitation of the two ions is also used 
for the luminosity measurement at RHIC |H1 E] • 

As the average impact parameter in such collisions is only about jTOj 
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one may expect strong field effects to be present for the pairs. 

Due to tfie design of tlie STAR detector only electrons and positrons having a trans- 
verse momentum pt > 65 MeV/c and being emitted with a rapidity \y\ < 1.15 can be 
detected. As most of the pairs produced in ultraperipheral collisions (UPCs) are emitted 
with energies of the order of a few mgC^ and at small angles, such a measurement can 
only look at the tails of the distribution of the pairs. Even this momentum range was 
only possible due to a lowered magnetic field in the STAR detector. 

One approach to the calculation of this process is the use of the impact parameter 
dependent equivalent photon approximation ^T1|^^1I2]- The additional electromag- 
netic processes are easily incorporated in this semiclassical approach and the cross section 
can be expressed as 

d^p^d^p. ~ J (Phduj^du^ d^p^d^p. • 

where cr^^^e+.e- denotes the cross section for real photons and d^L/ d%dujiduj2 the impact 
parameter dependent photon-photon luminosity (for the details of the photon-photon 
luminosity, see Sec. 2.7 of [T^ and references therein). One difficulty in this approach 
is the correct choice of the cutoff parameter present in the expression for the photon 
spectrum. This is especially difficult due to the smallness of me, which is much smaller 
than the "usual cutoff" imposed on the maximal transverse momentum of the photon 
from the elastic form factor of the ion, which is given by 1/ Ra ~ 80 MeV. It was found 
that the total cross section for electron-positron pair production is only reproduced with 
a cutoff chosen around uif.. On the other hand it was also found that neither a choice of 
TTLe nor of Ra is able to predict the total probability at impact parameter smaller than the 
Compton wave length Ac = 386 fm p. For a discussion about the choice of the impact 
parameter see, e.g., jTHj. In addition in the usual semiclassical approach the transverse 
momentum distribution of the photons is integrated over. In order to get, for example, 
the transverse momentum distribution of the pair, one needs to take this momentum 
distribution into account to get the correct final result. For a possible approach taking 
this into account from first principles see pT| IT^ . 

In Section El we show how our calculation is done in lowest order QED and in the 
semiclassical approximation. This is then used in Sec. El to calculate total cross sections, 
differential distributions for RHIC and also for possible LHC conditions. The comparison 
with the experimental results has been done in the meantime and will be presented by 
the STAR collaboration in another publication |19j . 

2 Calculation of pair production and nuclear excita- 
tion in lowest order QED 

The STAR experiment at RHIC measures the pair production cross section together with 
the double electromagnetic excitation in both ions, see Fig. ^ In order to incorporate the 
experimental conditions in the theoretical calculation, it is most appropriate to work in 
the semiclassical approach. Using the fact, that in this approach the probabilities of the 
individual processes factorize and are given by the product of the individual probabilities 
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Figure 1: The pair production together with the electromagnetic excitation of both ions, 
predominantly to the giant dipole resonance (GDR) is shown as one typical Feynman 
diagram. The process in lowest order involves at least the exchange of four photons 
(many more "soft Coulomb photons" are exchanged as well). Due to this the process 
predominantly occurs at small impact parameter, where the electromagnetic fields are 
strong, in contrast to the unrestricted cross section (without electromagnetic excitation 
of the ions), which has contributions coming also from large impact parameters. 



(for a theoretical description of this approach, see [^) we can write the cross section for 
this process as 

- 27r / hdhP^j^^ih)— — -— . (3) 



The minimum impact parameter was chosen to be hmin = '^Ra ~ 14fm, where we as- 
sume the nuclei to touch, that is, interact hadronically with each other. The use of the 
semiclassical description in this case is not only justified due to the strong Coulomb in- 
teraction between the two ions (leading to a large number of "soft" Coulomb photons 
exchanged between them), but also simplifies the calculation of this higher order process 
(pair production plus two GDR excitation processes) considerably. 

Following j2ilj the probability for GDR excitation in one ion is to a good approximation 
given as 

PcDRib) = S/b' (4) 

with 

S = " ^ 5.45 X lO-'Z'NA-^/Hm'. (5) 

where ttiat denotes the nucleon mass, and the neutron-, proton-, and mass-number of the 
ions are N, Z, and A respectively (we consider only symmetric collisions here, the calcu- 
lation can trivially be extended to incorporate also asymmetric systems). The excitation 
probability is inversely proportional to the energy lu (^ 80MeVA-i/3) of the GDR state. 
Neutrons are not only emitted from the GDR excitation but are also coming from higher 
excited states [S1I^I221- These can be taken into account approximately by increasing S 
accordingly. As this does not change the 1/6^ behavior for the small impact parameter, 
we are interested in, and only leads to a rescaling of the total cross section, not the form 
of the differential distributions, we have used in our calculation the more simple value of 
S in Eq. Of course one needs then to include in addition the different decay chan- 
nels into one, two, etc. neutrons j23l21ll2S], as well as, the fact that the probability for 
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GDR excitation is already quite large (about 0.5 for bmin) so that multiphoton excitation 
mechanisms need to be included as well. Assuming a Poisson distribution for the different 
(independent) excitation processes, one would need to replace Pgdr{P) then with [21] 

P{h) = l-eMPGDR{h)). (6) 

One sees that the multiphoton excitation tends to reduce again the probability (it has 
to stay below the unitarity limit of 1). In the appendix we show how a more complex h 
dependence than a simple 1/6^ for PcDRib) can be calculated within the same approach. 
Therefore an improved calculation can be done in this way. On the other hand for the 
current experimental accuracy the simplified approach seems to be appropriate. 

With this the cross section for pair production together with the GDR excitation of 
both ions is given by 

where we have introduced the two-dimensional Fourier transform of the impact parameter 
dependent probability for pair production d^P{h) / d^p+d^p- as 

We rewrite the integral over h in dimensionless units as 

I ^Jo{qb) =q^ f ^Mx) =: q"^ h{qbmin) ■ (H) 

Following the derivation of ^ |2E]; one can calculate the two-dimensional Fourier 

transform ^Sp^s'p in lowest order QED. The two Feynman diagrams for this process in 
the semiclassical approximation are shown in Fig. El One gets the differential probability 

as 

xTr|(^_ + m) [a^2D^ t/^)(^_- + m) ^/^^ + N^^ ip\4i- ^+ + m) 

X (^+ - m) [n^j^ 4[ + m) t/i) + N^^ ^^\4[- A + m) t/^)] |. (12) 

with 

N, = -ql N, = -{qf-{p++p^))\ 
N3 = -{qi + q)\ = -[gi + (g - p+ - p_)]\ 
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Figure 2: The two Feynman diagrams contributing to pair production in lowest order 
QED are shown. The crossed denote the couphng to the (external) Coulomb field of one 
of the ions. 

N2D = -{qi - P-f + m^, N2X = -{qi - P+? + 

N^D = -[qi + {q- p-)f + = -[qi + {q- p+)? + 

with the longitudinal component of qi given by gio = \ [(e+ + e_) + P{p+z + P-z)], qiz = 
^ [(e+ + e_) + P{p+z + P~z)] = ^qio and u^^'"^^ = 7(1,0,0, ±/?) are the four velocity of 
the two ions, 7 and f3 the Lorentzfactor and velocity of each ion in the center of mass 
frame. We have included in addition a nuclear form factor F{q). We choose for ease of 
computation in our case a monopole form factor of the form 

where = ^^2> is set to about 80 MeV in order to reproduce the rms radius of the 
ion. This leads in the terms Nq^, N{'^, Nj^'^ and N^'^ to a replacement of the term 1/g^ 
by F{q)/q^. The integration over (Pqi can be done analytically, using the usual tricks for 
Feynman integrations in two dimensions. For details of this, we refer the reader to ^J. 

The integral I^{z = qhmm) in Eq. (fTT|) can be solved analytically and calculated easily, 
as is shown in the appendix. 

Finally we make the integral over dq dimensionless to get 



d^^e+e-,2GDR ^ (2vr)2^ /" ^3^^ ( ) j. 



d^p+d^p- bt,inJ ~ ^~ d^p+d^p- ^^^^ 

For the unrestricted differential cross section (that is without triggering on the ad- 
ditional electromagnetic excitations of the ions), we integrate over h without the factor 
Pgdr{^) i^ Eq. (jni)- As the contribution coming from h < h„iin is small, we have extended 
the integration over all h. 

d'a,+,- f^,^ d^P{b) 



d^p^d^p^ J d^pj^d^p 



d d q — — expfzgO) 

d-^pj^d-^p_ 

.2 f,2^. d'Hq) 



{2'ay d'q 6{q} 



d^p+d^p_ 
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(15) 



This approach was pursued in [26 and total and differential cross sections were calculated. 
For the unrestricted cross sections we do not take a nuclear form factor of the two ions 
into account, as pair productions occurs predominantly at large impact parameter and 
for small of the two photons. 

Whereas for the unrestricted cross section only the value of P{q) for g = is needed, 
here our expression is a folding over a range of q given in terms of 1/bmin- In order to 
compare the differential distributions in both the restricted and in the unrestricted case 
we have also made calculations of the unrestricted cross section with the same kinematical 
restrictions as in the case with GDR excitations. 

The expression of Eq. ()14j) is in a form, which can be evaluated using a Monte Carlo 
integration for both the integration over z, as well as, the six-dimensional integration 
over and p_ at the same time. For this we have used VEGAS [21^. Both P{q) and Jq"* 
are oscillatory functions, having both positive and negative values, which could lead to 
cancellations. Looking at z^I{z) together with the result of P{z/hmin) one sees that the 
integrand falls off for large z, that is for large q. It is found that the main contribution 
comes from the region around z = 2 and that the contribution from the negative part at 
larger z are suppressed. The integration will have positive and negative contributions but 
the cancellations between them are not severe. With the help of VEGAS we can get the 
total cross section and also differential cross sections by binning the differential results. 

3 Results 

We have made calculations of the total cross section and differential distributions of the 
electron, the positron and the pair including the experimental restrictions at STAR. The 
integration over h (or equivalently q) is incorporated into the Monte Carlo integration. 
Another strategy would be to calculate d^P{q)/(Pp+d^p- for different values of q and 
fixed values of p+,p~ and do a Bessel transform in each case. For the total cross section, 
that is P{q), this can be done and was done as an independent check of our approach. 
To obtain differential cross sections, this approach is rather cumbersome. 

In a first step, we have calculated the total cross section as a function of bmin, the 
minimum impact parameter, by using three different approaches: We can calculate P{b) 
directly for the pair production process and integrate numerically over b. Second we 
can start from P{q) directly and do the integration over q numerically via the Fourier 
transformed of PcDRib), that is, using I{z). Finally we have done the calculation with the 
integration over q, that is z, directly with the Monte Carlo integration. In all three cases 
we have restricted the phase-space integration over the momenta of electron and positron 
according to the experimental conditions of STAR: pt > 60 MeV/c and \y\ < 1.15 for 
each lepton. The results are shown in Fig El The lines correspond to calculations with 
and without a monopole form factor for the nucleus, showing that the incorporation of a 
form factor is important. All three results agree quite well with each other, showing that 
our approach is working well. 

From this we get a total cross section for AuAu collisions at RHIC, including the 
restrictions \pt\ > 60 MeV/c, \y\ < 1.15 of 2.30, 1.76, 1.43 mb, for = 13, 14 and 
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Figure 3: The total cross section including kinematical restrictions of the STAR experi- 
ment are shown as a function of the minimal allowed impact parameter &mm- Calculations 
with and without a form factor for the nucleus are shown. The lines were calculated by 
calculating P{b) first and integrating then over b. These lines are in perfect agreement 
with a similar approach, where P{q) is calculated first and then integrated over q. The 
circles correspond to the results of our Monte Carlo approach, where, as explained in the 
text, both the integration over q and over p+ and p_ are done within the Monte Carlo 
integration routine. 

15fm, respectively. In addition we have calculated a number of differential distributions, 
which were also studied at STAR. The transverse momenta and energy distribution of the 
electron and positron are shown in Fig. |3 In lowest order QED the distribution of electron 
and positron are identical to each other. The difference between the two distributions can 
therefore be seen as a measure of the accuracy the MC integration. One can look also 
at properties of the produced pair: The transverse momentum and the invariant mass 
of the pair are shown in Fig. El We have not shown the rapidity distribution, which we 
found to be more or less flat over the allowed range. In all four diagrams we show also 
the differential distributions of the unrestricted cross section. We have rescaled the data, 
so that the total cross section are the same in both cases. We see that the transverse 
momentum distribution and the energy distribution of the individual leptons are more 
or less identical in shape. The same is also true for the invariant mass distribution, with 
the only exception that the "tail" at low invariant masses is higher for the unrestricted 
distributions. The biggest effect is seen in the transverse momentum distribution of the 
pair. Here we also expect the effect of the small impact parameter (corresponding to larger 
transverse momenta of the photons) to be largest. The slower fall-off at larger transverse 
momenta is most probably due to the fact that no form factor was used in the calculation 
of this cross section, which should be visible at Pt > 80 MeV/c. In order to investigate the 
effect of the nuclear form factor and in order to understand the large difference between 
the cross section with and without nuclear form factor, see Fig. 01 above, we show in this 
plot also the transverse momentum distribution of the pair for a calculation without form 
factor. It can be seen that in this case the cross section gets sizeable contributions for 
Pt > 80 MeV/c. 
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Figure 4: Energy and transverse momentum distribution are shown for the electron and 
the positron. In lowest order QED the two distributions are identical and the spread be- 
tween the two is a measure of the uncertainties coming from the Monte Carlo integration. 
This is compared with a unrestricted cross section calculation (dotted line), normalized 
to give the same total cross section. 



We have studied in addition the question, whether especially the distribution as a 
function of Pj of the pair is sensitive to the form of PcDRip) used in our approach. Using 
a more general approach, see appendix for details, one expect that the next correction is of 
the form Pgdr^^) ~ instead of PcoRi^) ~ 1/^^- such a model the average impact 
parameter changes from about ^ 19 fm to about ^ 17 fm, which is a small 
change compared to the Compton wavelength of the electron (400 fm), but is still a 10% 
reduction of this average impact parameter. We have therefore studied the distribution 
of electron and positron for this restriction. As expected the differential distributions are 
found to be the same within the uncertainties of the Monte Carlo approach used. 

As an outlook for future experiments we are showing results and distributions one 
might expect to see for PbPb collisions at the LHC. Using the same kinematical restric- 
tions as for the STAR experiment, the results are shown in Fig. IHl and [3 

As a rather optimistic estimate we have calculated also the differential cross section 
for a kinematical range of \pt\ > 2.6 MeV/c and \y\ < 1.5, where ALICE will be able to 
detect the electrons with its Internal Trigger System (ITS), even though it will not be 
able to measure energies or momenta. 

With these kinematical conditions we can study the question, whether ALICE will be 
able to see multiple pairs produced in a single collisions. For this we calculate the impact 
parameter dependent cross section under the kinematical conditions. One finds that for 
impact parameters close to bmin Pip) ~ 20%. Following [2H1 EHl EHl EDI we use a Poisson 
distribution as a good approximation for the pair production probability 

m,&) = -y-exp(-P(6)) (16) 

Multiplying with PcD^ib) and integrating over b we get the cross section for one, two, 
. . . pair production. The result as a function of bmin is shown in Fig. |H1 together with the 
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Figure 5: The invariant mass and the transverse momentum of the e~^e~ pair are shown. 
Again the results are compared with the unrestricted cross section (dotted hne). In addi- 
tion we have plotted the result of the restricted calculation without form factor (dashed 
line). Both results have been normalized to give the same total cross section. 



total cross section 

oo 

O'total = 

E (17) 

N=l 

The cross section one would get from the Born cross section can be interpreted as a 
"multiplicity" cross section 



/oo 
d%P\GDR, b)P{b) = J2 <^iN) (18) 
AT 1 



N=l 

and would be relevant in order to calculate the number of pairs produced (in contrast 
to the number of events). One can see that about 10% of all events are multiple pair 
production events and accordingly also about 10% of all pairs are produced in a multiple 
pair production process. This shows that at ALICE one should be able to detect and 
study multiple pair production. 

We have investigated a similar question also for RHIC using as an estimate for a 
possible range pt > 50 MeV/c and 2.5 < y < 4.0. Unfortunately the probability for pair 
production under these conditions is only of the order of a few permille and therefore 
the multiple pair production cross section is less than one permille of the single pair 
production cross section, making such an investigation difficult. 



4 Discussion and Outlook 

We have calculated total and differential cross section of the pair production process 
in ultraperipheral heavy ion collisions in lowest order QED for the simultaneous elec- 
tromagnetic excitation of both ions. We have seen that the most sensitive quantity is 
the transverse momentum distribution, which differs mostly from the distribution of the 
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Figure 6: Energy and transverse momentum distribution are shown for the two leptons 
for PbPb colhsions at the LHC (sohd hne). This is compared with the rescaled spectrum 
of AuAu colhsions at RHIC (dotted hne). 



unrestricted cross section. As the comparison wiU show [13 12] our results were found 
to be in good agreement with the experimental results. On the other hand only about 
50 events were found at STAR, so the overall statistics is not very good. Additional 
runs might give better statistics. Still our analysis shows that the data at the moment 
give no sign that higher order Coulomb effects are large for pairs produced with these 
large transverse momenta [21] . Such higher order Coulomb effects would lead most likely 
to an asymmetry of the electron and positron distribution especially for the transverse 
momentum distribution. At the moment however no calculation exists, which describes 
consistently the effects of the strong Coulomb fields of both pairs on the pair production 
process at small impact parameter. The experimental condition of RHIC does not allow 
to look for multiple pair production effects as the probabihty for pair production with 
these conditions is well below one. 

As already mentioned above the transverse momentum cut of pt > 65 MeV/c at STAR 
was only possible due to a reduced magnetic field. There are currently plans to use even 
lower magnetic fields and also making use of other detectors within STAR [32 in order 
to extend the measurements both to smaller transverse momenta and to larger rapidities. 
It remains to be seen, whether the new phenomenon of multiple pair production, will 
then be detected. On the other hand with the low transverse momentum cutoff of the 
ITS at ALICE/LHC, about 10% of all pair production events are going to be multiple 
pair production processes, therefore one should expect that this new phenomena will 
be observed easily there. The fact that no kinematical information and also no particle 
identification is possible at ALICE, will make such a measurement still a challenge. 
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The invariant mass and the transverse momentum of the e^e~ pair is shown 
colhsions at the LHC (sohd hne). The resuhs are compared with the rescaled 
of AuAu colhsions at RHIC (dotted hne). 
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in understanding the possibihties of detecting the pair production at LHC. 



6 Appendix 

Throughout our calculations we have assumed that only the GDR excitation is relevant 
for the triggering and therefore a simplified dependence on 5, see Eq. (jl)) has been used. 
In this appendix we want to show that this is not a real limitation, but that other 
impact parameter dependencies can be treated as well. For example assuming that the 
higher resonant states of the GDR are excited through a Poisson process, we would 
need to replace PcDRip) by 1 — exp(— Pg'£)/j(6)). In general we assume that the relevant 
PA^A*^x+xnip) can be expressed as a series of inverse powers of h 

OO Q 

PA^A'^X+Uh) = T.l^ (19) 
n=0 " 

Using this in the expression for (fae+e- ,2Gdb./ (d^P+d^P-) , see Eq. (|T^ . we need to calcu- 
late generalizations of I{z), Eq. (Illj) of the form 

^"(^) ■= r s^o(^) (20) 

with hi^z) corresponding to the one used in our calculations. 

For the calculation of these integrals, we first use the following recursion relation 

r(,^-_lM -IM (,] (o^] 

^ ~ (n- 1)^-1 (n - l)2^"-2 (n - l)2'"-2^^^ 

This relation can be easily derived by partial integration and by using the well-known 
recursion relations between the Bessel functions Jn{z) By repeated application of 
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Figure 8: The different cross sections for single and multiple pair production are shown 
for the kinematical conditions at ALICE. The cross section is shown as a function of 
bmin, the minimal impact parameter. See the text for further details of the different cross 
sections. 



Eq. every with odd and even n can be reduced to the starting values Ii{z) and 
Io{z), respectively, which are given in the literature 



Jo(z) = l-;^iFji;l,|-^) (22) 



and 



h{z) = J2Fs\ 1, 1; 2, 2, 2; -- - In - - 7 (23) 



where 7 is the Euler constant. 

These two expressions can be calculated easily by the rapidly converging power series 
of the hypergeometric functions or by using suited polynomial expressions given in the 
literature jSSj- We are here only interested in the case, where n = 2m + 1 is an odd 
number, in which case the complete recursion relation is given by 

(_-[)m ( m , A.s r ™ / A \^ \ 

WW = ;;y;^{^o(.)E3!(^-l)!(-^) -2^iWE(^-l)!^(-;^) 

(24) 

Moreover this equation can be further simplified, by splitting off the terms singular at 
z = (the principal part of the Laurent expansion in z\ by using the power series for 
the Bessel functions: 



m 00 (_^2l^Y 



^..w-y £s(iTi^ (25) 

and by rearranging the resulting sum. After some straightforward algebra we thus obtain 
the compact expression for m > 

hm^iiz) = 2, J+ i)!2 ^' ^-^3 [l, 1; 2, 2 + m, 2 + m; --j (26) 

' % ^ + 7 ) + E (27) 



22™(m)!2 V 2 



A;=0 
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where 

«™.o = ^^^E-, = ^^(^^^z^y:;:!^;^^^^ A: > 1 (28) 
For m = 1 we get the exphcit expression 

Hz) = ^ + i (l^ i + 7 - l) - Y^^' 2^3(1, 1; 2, 3, 3; -z^A) (29) 

and similarly for m = 2. Again the hypergeometric function in Eq. ()27j) can well be 
calculated numerically by its power series: 

(l. 1; 2. 2 + m. 2 + -^j = (,„ + 1)P (-1) | (30) 

We use this rapidly converging power series in our numerical calculations. For the case 

where n is an even number, the same approach can be used. For completeness we only 
give here the final result: 

. / X 1 ^ f 1 3 (_i)"^22'"m!2 ^ ^ 

hmiz) = -7; 7z r 1F2 -m + -: -m + -, 1: + ^ 7 — . (31) 

; 22™-i(2m-l) I 2' 2' ' 4; (2m)!2 ^ ' 
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